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Abstract 

Motivated by a two-bump (or 1-peak plus 1-hump) structure in the ATIC data, we perform a 
statistical analysis fitting the PAMELA and ATIC data to a dark matter model, in which the dark 
matter particle can undergo both annihilation and decay. Using a chi-square analysis we show that 
both data can be simultaneously fitted better with such a double-action dark matter particle. We 
use an existing neutrino mass model in literature to illustrate the idea. 
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I. INTRODUCTION 



The year 2008 had been filled with excitement from a number of dark matter (DM) 
experiments. The PAMELA Collaboration [l| has reported an unexpected rise of positron 
fraction at the energy range of 10 — 100 GeV, unlike the power-law falling background. 
This provides further support to the earlier results reported from HEAT 2| and AMS- 
1 [sl. However, similar enhancement of the anti-proton fiux was expected but not seen 
by PAMELA |^ provides a challenging puzzle. Other surprises came from two balloon 
experiments ATIC js] and PPB-BETS 6| at the South Pole Antarctica. The ATIC data 
showed an excess of galactic cosmic-ray electrons/positrons at energies of 300 — 800 GeV. 
These experimental results have stimulated a lot of theoretical speculations about possible 
mechanisms, including dark matter annihilation [3], decaying dark matter with a very long 
lifetime js], or simply astrophysical origins from either ultrahigh energy cosmic rays 9] 



or nearby pulsars 



101 ] within a few kilo-parsec. If the observed positron excess is indeed 



due to dark matter annihilation, the data sets require an annihilation cross section {av) of 
the order of 10"^'^ cm^ s~^, which is two to three order of magnitudes larger than naively 
expected from a thermo-WIMP dark matter in most popular models like the minimal SUSY 
and Kaluza-Klein models. Either a large boost factor or Sommerfeld-type enhancement 11 1 
can be used to explain such a large annihilation cross section. On the other hand, a very 
long lifetime of the decaying dark matter of the order of 10^^ seconds is required to fit the 
data. Such a long-lived dark matter is consistent with other cosmological constraints on 
our Universe. Implications for further investigations in future gamma-ray experiments and 



neutrino telescopes have been studied 



12|. 



We point out that the excess in the ATIC data in fact consists of 1-peak plus 1-hump 
structure. The peak is from 300 — 800 GeV while the bump from 80 — 300 GeV. A possible 
explanation is that the dark matter particle in the Universe undergoes both annihilation 
and decay. The annihilation gives rise to the peak around 600 GeV while the decay is 
responsible for the small bump around 80—300 GeV. The resulting positron fraction observed 
at PAMELA is a combination of annihilation and decaying contributions. 

In this work, we perform a analysis which shows that the fitting using a single mode 
(either annihilation or decay) is far less satisfactory than the double mode (both annihila- 
tion and decay). The study shows that a dark matter particle of about 640 GeV with a 



2 



monochromatic annihilation spectrum and a soft decaying spectrum is the best simultaneous 
fit to the PAMELA and ATIC data. 

The PAMELA data is highly restrictive on the anti-proton mode of the dark matter 
annihilation or decay {4]. It points to the hint that the dark matter particle may be lepto- 
philic or carrying a lepton number. There are some models that propose the TeV right- 
handed neutrino, which is responsible for neutrino mass, to be the dark matter candidate. 
We will borrow an example in literature 
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15| to illustrate the possibility. 



Some of the highlights in this paper include 

• We show in a more quantitative way how the models are fitted to the ATIC and 
PAMELA data simultaneously; 

• We use the MINUIT program from the CERN library to vary the parameters of the 
model to minimize the in order to obtain the best fit of the model; 

• We show that the dark matter candidate that can both annihilate and decay can fit 
better to the data. 

The organization of the paper is as follows. In the next section, we summarize the 
formulae we employed in our analysis for the annihilation and decay of the dark matter 
particle without referring to any particular dark matter models. In Sec. Ill, we show our 
numerical analysis of fits. In Sec. IV, we describe a slight modification of an existing 
dark matter model in literature that may give rise to the double-action. We conclude in 
Sec. V. 



II. DOUBLE-ACTION DARK MATTER 



A. Dark Matter Annihilation 



Assuming a steady state condition while solving the diffusion equation for the positron 
as it traversed across the universe, its flux at Earth can be casted into the following semi- 

<^eAE) = '^feAE), (1) 



analytical form [16|, ll 
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with Ve+ close to the velocity of light c and the function /e+ (E) is given by 

1 rE„- 
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with -Emax = ^dm in the case of annihilating DM, Mdm and pdm are the mass and the density 
of the dark matter respectively, and the overall constant B is the boost factor. In Eq.([2]) we 
have expressed the source term Qann according to 



M, 



dm ' 



dE,_ 



(3) 



where t] = 1/2(1/4) for Majorana or Dirac particle. The summation is over all possible 
channels that can produce positron in the final state, and dN(.+ / dEf,+ denotes the spectrum 
of the positron energy per annihilation in that particular channel. We have suppressed the 
index labeling the various channels to avoid notation cluttering in the equations. The halo 
function I{\d) can be parametrized by 



I{\d) = Oq + Oi tanh 



'b,-V 



a2 exp( 



(l - b2f 



+ ^3 



(4) 



with / = log^Q (Az)/kpc) and the diffusion length Xd{E, E') is given by 

'{E'/QeVy-^ - (E/GeVy-^' 



d-1 



(5) 



The constants 00,1,2,3? &i,2; ci,2 and 6, Kq can be found respectively in Table 2 and Eq.(ll) of 



Ref. 



. The energy loss rate function b{E) in Eq.([2]) is 

b{E) 



(6) 



(GeV X te) 

where te = 10^^ seconds. In our analysis, we simply employ the monochromatic annihilation 
spectrum 



dN^ 



-{E) 



1 



-5 1 



E 



M, 



dm 



(7) 



Analogous formulas for the electron will be omitted here. We only consider the monochro- 
matic e~e^ spectrum in this study for simplicity and clarity. We could also add the muon and 
tau channels in the annihilation, but that would introduce more parameters to complicate 
our analysis. 
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B. Decaying Dark Matter 



The source term for a decaying dark matter in a particular channel is 

1 / Pdm \ dN,+ 

Wdec - — I , loj 

Tdm VMdm/ dhe+ 

where Tdm is the lifetime of the DM and dNe+ / dEf,+ is the positron energy spectrum per 
decay of the DM. The function ff,+ {E) is now given by 

fe^E) = ^ dE'I{Xo{E,E'))Q^UE') 

with -Emax = for the decaying dark matter and summation over all decay channels is 

explicit. The flux is the same as in Eq.([l]). 

In our analysis, we use either a (i) monochromatic decaying spectrum: 

^'^'*^E)^^s(l-m, (10) 



dEe+ A/dm V Mdm 

or (ii) varying decaying spectrum 

dNe. 80E / 2EY .... 

The exact form of the varying decaying spectrum is not crucial in our analysis. As long as 
it is soft it suffices to suit our purpose. We will show that the model that we will use in 
this work gives an energy spectrum consistent with the above varying spectrum. Analogous 
formulas for the electron will be omitted as before. 



C. Background fluxes 

The background electron/positron fluxes from astrophysical sources are believed to be 
mainly due to Supernova explosions for the primary electrons and from the interactions 
between the cosmic ray nuclei (mainly proton and helium) and atoms (mainly hydrogen 
and helium) in the interstellar medium for the secondary electrons and positrons. They are 



commonly parametrized as 



19| 



4.5E' 



0.7 



<T)bkdg ^ — /-|2^ 
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1 + 11^0-9 + 3.2^2-15 1 + 110^1-5 + 580^4-2 
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where E is in unit of GeV and the unit for the flux is GeV~^ cm~^ s~^ sr~^. We use a 
normahzation of 0.7 so that the flux calculation is consistent with the ATIC data in the low 
energy range of 20 — 70 GeV. ^ 



D. Propagation models 

To evaluate the halo function /(A/5), we will use the Navarro- Frenk- White (NFW) ha 
profile with the propagation models "M2", "Med" and "Ml" as parametrized in Ref. 18 | 
The models "M2", "Med" and "Ml" are characterized by the propagation length in the 
increasing order. The analysis we will perform in the next section can be straightforwardly 
repeated for the other popular halo models specified by the Moore profile and the cored 
isothermal profile. 



III. ANALYSIS 



Theoretical predictions of the energy spectrum and the positron fraction depend in general 
on the three input parameters Tdm, {cv) and M^m- To achieve the output more or less 
consistent with the ATIC and PAMELA data, we use Tdm = 1-3 x 10^'^ s, B{av) = 5.4 x 
10~^^ cm^ s~^ and = 643 GeV. The energy spectrum for ATIC and the positron fraction 
for PAMELA are shown in the Fig. [1] and Fig. [2] respectively. We have used the NFW halo 
model with propagation model "Med" in the figures. From these two plots it is interesting 
to see that both PAMELA and ATIC data can be fitted simultaneously by the double-action 
dark matter. We now want to justify this fact more quantitatively using the technique of 

fits. 



A. Fitting with PAMELA data only 

The PAMELA data shown in Fig. [2] has a rising trend starting from point #9 to #16 
which we will be using exclusively in our numerical analysis. The first eight data points are 



The normalization of background can vary between 0.6 and 0.8 so that our conclusion of 1-peak-l-bump 
structure is still valid. If we also vary the normalization in the fits, we found that the best fit for the 
ATIC data (M2) is similar to the result in row #4 in Table II and the normalization is 0.69, which is close 
enough to our fixed value of 0.7. 
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FIG. 1: The spectrum for the ATIC data. The dark matter mass is taken to be 643 GeV, with 
a monochromatic spectrum for DM DM — > e~^e~ annihilation and a soft decaying spectrum of 
dN/dE ~ E{1 - 2£;/Mdm)V^dm for DM ^ e+e"X decay. 
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■ magnetic polarity state 
l|. Both horizontal and 
but only the latter one 



of low energy less than about 10 GeV, where the solar activity o: 
is expected to play a significant role in the positron abundance^ 
vertical errors are explicitly given in the PAMELA publication 
will be taken into account in our analysis. The parameters used in the analysis are Tdm, {(tv), 
and Mdm as mentioned earlier. As indicated in previous section, we use a monochromatic 
electron/positron spectrum for dark matter annihilation, while for decaying dark matter 
we use either a (i) monochromatic (mono) or (ii) varying spectrum (var). It is denoted 
by "mono" or "var" in the tables. However, we emphasis that the results do not depend 
strongly on the exact spectrum as long as it is soft enough. We will show that in some cases 
the soft spectrum actually fits better than the monochromatic spectrum. 

In Table m we show the fits to the PAMELA data (point #9 to ^^16) with one or two of 
the parameters fixed. It can be seen that the "Ml" propagation model fits slightly better 
than "Med" but much better than "M2" . When we vary the mass of dark matter from 200 
GeV to 1000 GeV, the goodness of the fit, measured by the x^y is roughly independent of 
the mass. (It shows a slight better fit when increases, but not of any significance.) 
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FIG. 2: The positron fraction for the PAMELA data. The dark matter mass is taken to be 
643 GeV, with a monochromatic spectrum for DM DM — > e^e~ annihilation and a soft decaying 
spectrum of dN/dE ~ E{1 - 2E/Mdmf/Ml^ for DM e+e~X decay. 

Since the PAMELA data did not show any peak structure, the data would not prefer any 
mass of the dark matter. As long as the propagation diffuses the spectrum or the spectrum 
itself is soft, it can fit the data well. It is shown in the first three rows in Tabled when we 
used the varying spectrum. The "var" spectrum fits much better than the "mono" spectrum. 
Thus, the PAMELA data alone do not constrain the mass of dark matter to any significant 
extent according to this analysis, as long as the mass is heavier than about 200 GeV. 



B. Fitting with ATIC data only 

The ATIC data exhibits a more interesting feature of one peak plus one hump structure. 
Thus, we expect that the data prefer some mass range of the dark matter for both annihi- 
lation and decaying contributions. By comparing row #4 with #1 to #3 in Table HTl where 
the "M2" propagation model is used, we note that with a dark matter about 744 GeV the fit 
is substantially better with both annihilation and soft decaying contributions included than 
just either one of them is used. On the other hand, the other two propagation models do 
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not show such effects. The less diffuse propagation model "M2" fits better than the "Med" 
and "Ml" models. 

C. Fitting using both PAMELA and ATIC data 

The fits are shown in Table IIIII Similar to the fits with PAMELA data only, the fits 
using "Ml" propagation model are slightly better than using "Med", which are in turns 
much better than "M2". The fits using "M2" is not good at all, given the fact that per 
d.o.f. is large. This behavior is similar to those fitted to PAMELA data alone. 

In the fits with "Med" propagation model, we first fitted with a negligible decaying 
contribution, i.e., with annihilation contribution only. The best that we can do is — 59 

TABLE I: Fitting with PAMELA data only. When we fix Tdm = 10^° s, the contribution from 
decaying is negligible. Similarly, when we fix B{av) = 10"^'^ cm^ s^^ the contribution from 
annihilation is negligible. The "mono" means monocliromatic spectrum, whereas "var" means 
varying spectrum for the decaying contribution. 



Prop. 


''"dm 


B{av 


) 




XV# d.o.f. 


Comments 


model 




(cm^ 




(GeV) 








10^0 (fixed) 


10-23 




300 (fixed) 


79.6/7 


mono 


NEW M2 


0.45 • lO^^^ 


10-40 


(fixed) 


212 


18.1/6 


var 




0.55 • lO^^^ 


0.77- 


10-24 


200 (fixed) 


10.1/6 


var 






0.15 • 


10-23 


250 (fixed) 


3.0/7 


mono 




fixed 


0.21 • 


10-23 


300 (fixed) 


2.9/7 


mono 


NEW Med 


at 


0.37- 


10-23 


400 (fixed) 


2.8/7 


mono 




10« 


0.83 • 


10-23 


600 (fixed) 


2.7/7 


mono 






0.15 • 


10-22 


800 (fixed) 


2.7/7 


mono 






0.22 • 


10-22 


1000 (fixed) 


2.6/7 


mono 




fixed 


0.59 • 


10-24 


200 (fixed) 


3.6/7 


mono 


NEW Ml 


at 


0.36 • 


10-23 


500 (fixed) 


2.2/7 


mono 




10« 


0.70 • 


10-23 


700 (fixed) 


2.0/7 


mono 






0.14 • 


10-22 


1000 (fixed) 


1.9/7 


mono 
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at Mdm — 600 — 700 GeV. When we also turn on the decaying contribution with a varying 
spectrum, the goes down to 42. This result supports a dark matter that annihilates and 
also decays with a soft varying spectrum. The best fit is 

Tdm = 0.13 ■ 10^8 s, B{av) = 0.54 • lO'^^ cm=^ s-\ M^m = 643 GeV (14) 

which has a — 42/26 d.o.f. We used this set of fitted parameters in Figs. [T]and[2l 

The fits using the "Ml" propagation model improve further. Features are similar to the 
case of "Med". The best fit using monochromatic annihilation and decaying contributions 
is (second last row of Table IIIip 

Tdm = 0.38 ■ 10^^ s, B{av) = 0.61 ■ 10^^^ cm^ s~\ M^m = 745 GeV (15) 

with a = 38.5/26 d.o.f. The best fit using monochromatic annihilation but varying 
spectrum for decaying contribution is (last row of Table IIIB 

Tdm = 0.33 ■ 10^9 s, B{av) = 0.54 ■ lO'^^^ cm^ s-\ M^^ = 643 GeV (16) 

with a = 38.6/26 d.o.f. 



TABLE II: Fitting with ATIC data only. Other details are the same as Table HI 



Prop. 


Tdm 




B{av) 




XV# d^o.f. 


Comments 


model 


(S) 




(cm^ s""*^) 


(GeV) 








1040 


(fixed) 


0.51 • 10-23 


536 


39.5/19 


mono 


NFW M2 


0.36 


•1027 


10-40 (fixed) 


1072 


39.5/19 


mono 




1.0- 


1026 


10-40 (fixed) 


3190 


25.3/19 


var 




0.65 


•1027 


0.80 • 10-23 


744 


23.7/18 


var 




1040 


(fixed) 


0.78 • 10-23 


745 


27.4/19 


mono 


NFW Med 


0.32 


•1027 


10-40 (fixed) 


1490 


27.4/19 


mono 




0.25 


•1029 


0.78 • 10-23 


745 


27.4/18 


var 




1040 


(fixed) 


0.68 • 10-23 


740 


34.2/19 


mono 


NFW Ml 


0.36 


•1027 


10-40 (fixed) 


1470 


34.2/19 


mono 




0.11 


•1027 


10-40 (fixed) 


4420 


37.6/19 


var 
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TABLE III: Fitting with PAMELA and ATIC data. Other details are the same as previous tables. 



Prop. 


7"dm 


B{av 


) 




XV# d-o.f. 


Comments 


model 


(s) 


(cm^ 




(GeV) 








10^0 (fixed) 


0.31 • 


10~23 


400 (fixed) 


297/28 


mono 




10^0 (fixed) 


0.47- 


10-23 


500 (fixed) 


281/28 


mono 




10^0 (fixed) 


0.70- 


10-23 


600 (fixed) 


269/28 


mono 


NEW M2 


10^0 (fixed) 


0.94 • 


10-23 


700 (fixed) 


275/28 


mono 




10^0 (fixed) 


0.12 • 


10-22 


800 (fixed) 


284/28 


mono 




10^0 (fixed) 


0.16- 


10-22 


1000 (fixed) 


343/28 


mono 




10^° (fixed) 


0.55 • 


10-23 


535 


267/27 


mono 




10^° (fixed) 


0.27- 


10-23 


400 (fixed) 


96.8/28 


mono 




10^° (fixed) 


0.42 • 


10-23 


500 (fixed) 


74.3/28 


mono 




10^° (fixed) 


0.60- 


10-23 


600 (fixed) 


59.1/28 


mono 




10^° (fixed) 


0.80- 


10-23 


700 (fixed) 


59.3/28 


mono 


NEW Med 


10^0 (fixed) 


0.10 • 


10-22 


800 (fixed) 


63.2/28 


mono 




10^° (fixed) 


0.15 • 


10-22 


1000 (fixed) 


102/28 


mono 




0.16 • 10^8 


0.62 • 


10-23 


745 


55.9/26 


mono 




0.13 • 10^8 


0.54- 


10-23 


643 


41.9/26 


var 




10^0 (fixed) 


0.34- 


10-23 


500 (fixed) 


55.3/28 


mono 




10^0 (fixed) 


0.48- 


10-23 


600 (fixed) 


41.1/28 


mono 




10^° (fixed) 


0.64- 


10-23 


700 (fixed) 


38.8/28 


mono 




10^° (fixed) 


0.83- 


10-23 


800 (fixed) 


40.0/28 


mono 


NEW Ml 


10^0 (fixed) 


0.99 • 


10-23 


900 (fixed) 


66.9/28 


mono 




10^0 (fixed) 


0.12 • 


10-22 


1000 (fixed) 


66.7/28 


mono 




0.38 • lO^s 


0.61 • 


10-23 


745 


38.5/26 


mono 




0.33 • 10^9 


0.54- 


10-23 


643 


38.6/26 


var 
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IV. A MODEL 



We use a TeV right-handed neutrino mass model to present a dark matter candidate 
that can annihilate and decay. The model can be described by the following interaction 
Lagrangian 

Ant = faf3LlCiT2LpS^ + giaNiS^iaR + g2aN2S:^iaR + 'ii-C. 

+ Mn.N^CNi + Mm,N^CN2 - V{Si, S2) (17) 

where and laR are the lepton doublet and singlet respectively with a, (3 denoting the 
family indices, A^i 2 are the two right-handed neutrinos, C is the charge-conjugation oper- 
ator, and V{Si, S2) is the scalar potential for the two complex scalar fields 5*1 and S2 that 
containing a term \s{Si S2Y ■ Note that fnp is antisymmetric under the interchange of the 



family indices. The original model in Ref. [13| contains only one A^^i, but the improvement 
in Ref. jl^ by adding another A^2 niakes the model consistent with the neutrino oscillation 
data. For simplicity we shall only describe the lighter right-handed neutrino A^^i (denoted 
simply by in what follows), which is in the TeV mass range, as far as dark matter is 



concerned. It has been shown ij] that can be a dark matter candidate. The annihilation 
NN —>■ e~^e~ can go through t- and w-channel diagrams with an intermediate 5*2" (see Fig. 
|3]^a)) with the rate given by 

CrVrel = — / dx < ——5 —5 — --r + 



647r s J -I [4{Ml - M| - f (1 - P^x))^ 4{Ml - - f (1 + (3p,x))^ 

'^_Mk 1 ri8) 

(Ml-Ml-l{l-P^x)){Ml-Ml-l{l + P^x))j ^ ' 

where (3^ = (1 — 4M^/s)^^^. As (3^ — 0, the above annihilation rate vanishes. This is 
expected for the annihilation rate for identical Majorana fermions is P-wave suppressed 



201 1 . When the center of mass energy ^/s is slightly above the threshold of 2Mm, the 
annihilation electron/positron energy spectrum is almost a monochromatic one. We show 
in Fig. m the annihilation cross section versus the right-handed neutrino mass for various 
values of the scalar mass, taking the center-of-mass energy y/s at which Vj-e\ ~ 10""'^. The 
plot demonstrates the P-wave suppression. Therefore, a large boost factor of order O(IO^) is 
needed to fit the data in this model. However, this suppression might not be taken literally. 
If one is willing to extend the model by introducing new long range force among the dark 
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matter, the Sommerfeld enhancement for P-wave annihilation can be significantly larger 
than the S-wave case 2l[ |. 




FIG. 3: Feynman diagrams for (a) annihilation NN ^ e and (b) decay N ^ e e'^Vf 
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FIG. 4: Annihilation cross section versus the right-handed neutrino mass for various values of the 
scalar mass at a center-of-mass energy y/s when the Urei ~ 10"'^. 

In Refs. [isl, [l^ , the is assumed stable by imposing a Z2 parity. Here we introduce a 
small violation of this parity by adding a term eSiS2 + H.c. to the scalar potential, where 
e ~ (1 eV)^, which is of the order of the square of neutrino mass. The decay of can then 
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go through a Feynman diagram shown in Fig. [3]^b). ^ An order of magnitude estimate of 
the decay width of the can be given as 

~ glJl^e^Mn/Mt^ . (19) 

Taking typical values of the couplings [l^ {gie ~ 10^^, /12 ~ 10"^, Mn ~ Ms^ ~ TeV), the 
lifetime of is roughly 

rj, ^ 1026 1^^^ sec . (20) 

It is interesting to see that when e is of the order of the neutrino mass squared, the amount 
of violation of the Z2 parity is in the right order to fit the data. We calculate the normalized 
energy spectrum of the decay A^ — > e~e~^i'^ shown in Fig. [5l where the approximate spectrum 
80x(l — 2xY with X = E/M^m is also shown. The figure justifies the approximation of the 
energy spectrum that we have used in our analysis given in the previous section. Note that 
the exact form is not crucial in the fits as long as the spectrum is soft. In this model, the 
decay and annihilation of the dark matter are pure leptonic. It will not give enhancement 
to the p flux. 



^ There is another decay channel N ~> e fi^Vf. given the same couphngs. The hfetime will be shorten by 
a factor of two, but it does not affect our ordcr-of-magnitude estimate. 
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0.1 0.2 0.3 0.4 0.5 



X = EM^^ 

FIG. 5: Normalized energy spectrum 1/TdT/dx with x = E/Mdm for the decay N — > e~e~^i'^ and 
the approximation 80x(l — 2x)^. 



V. CONCLUSION 

We have pointed out a 2-bump (or 1 peak plus 1 hump) structure in the ATIC data. If 
this feature can be sustained it cannot be explained by dark matter annihilation or decay 
alone. It can be either a two-component dark matter or the dark matter can undergo both 
annihilation and decay at the same time. We have shown that such a double-action dark 
matter can fit better to the ATIC and PAMELA data simultaneously than just annihilation 
or decay alone. We have employed a TeV right-handed neutrino model to illustrate the 
idea. The original model only allows dark matter annihilation, but here we have introduced 
a small breaking of the Z2 parity at the order of the neutrino mass squared. With such 
a coincident size of Z2 breaking one can explain the long lifetime required to fit the data. 
This indicates there may be intricate connection between neutrino mass problem with dark 
matter physics. We close with a few comments. 

• The model of TeV right-handed neutrino can account for the neutrino mass and oscil- 



lation pattern 
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14i\ . It is also consistent with lepton-flavor violation. The lightest 
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of the right-handed neutrino can be a dark matter candidate and its rehc density can 
account for the observed dark matter density. 

• Reference jl^ showed that the model can be made consistent with existing neutrino 
oscillation data by tuning the parameters of the model. If so there would also be the 
muon and tau channels in the annihilation. The electron/positron spectrum would 
contain the monochromatic part and the continuous part from muon and tau decays. 
The analysis would be much more involved and beyond the scope of this paper. 

• The decay and annihilation of the dark matter are pure leptonic. It will not give 
enhancement to the p flux. 

• The Z2 parity violation parameter e that we introduced is, by accident, at the order 
of the square of the neutrino mass. 

• The prediction for the gamma-ray flux mainly comes from the bremsstrahlung off the 
electron / positron. 

• Since the TeV right-handed neutrino has no appreciable coupling to quarks or gluons, 
the scattering cross section with nuclei is negligible. Thus, the sun or the Earth will 
not capture any large amount of dark matter, and so no enhancement to the neutrino 
flux coming from the core of the sun. 

• Either annihilation or decay alone cannot explain the probable 1-peak plus 1-hump 
structure in the ATIC data. We have shown that a dark matter that can annihilate 
and decay simultaneously can explain the ATIC and PAMELA data at the same time. 

Dark matter interpretation for the PAMELA and ATIC experiments is exciting since it 
implies new physics beyond the Standard Model. However, one should keep in mind that 
a yet unidentified astrophysical object such as a nearby pulsar or micro-quasar could be 
a primary source as well. Furthermore, many models are capable to explain the excess 
anomaly. PAMELA is extending the spectra measurement to higher energy of about 300 
GeV for positrons and 500 GeV for electrons. New data is also expected soon from the 
FERMI satellite for the diffuse Galactic cosmic 7-ray spectrum. These future developments 
will certainly help to discriminate models and unravel the true nature of the anomaly seen 
thus far. More excitements are waiting ahead of us! 
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